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We describe a tensorial generalization of the Navier slip boundary condition and illustrate 
its use in solving for flows around anisotropic textured surfaces. Tensorial slip can be de- 
rived from molecular or microstructural theories or simply postulated as an constitutive 
relation, subject to certain general constraints on the interfacial mobility. The power of 
the tensor formalism is to capture complicated effects of surface anisotropy, while pre- 
serving a simple fluid domain. This is demonstrated by exact solutions for laminar shear 
flow and pressure-driven flow between parallel plates of arbitrary and different textures. 
From such solutions, the effects of rotating a texture follow from simple matrix algebra. 
Our results may be useful to extracting local slip tensors from global measurements, 
such as the permeability of a textured channel or the force required to move a patterned 
surface, in experiments or simulations. 



1. Introduction 

The emergen ce of microfiuidics has focused renewed attention on hydrodynamic bound- 
ary conditions ( Stone et al. 2004 ). R educing fluid volumes enhances the impact of surface 
phenomena ( Squires fc Quakef 20051 ) . so the use of appropriate boundary conditions is 
crucial to the design and optimization of lab-on-a-chip devices. It is now widely recog- 
nized that the classical no-slip hypothesis supported by macroscopic experiments does 
not always apply at the micro- and, especially, the nano-scale. 

In this context, the phenomenon of liquid slip at solid surfaces h as been studied ex- 
tensively in experiments, theoretical calcu lations, and simulations ( Vinogradova 19991 : 
Lauga et aZ.ll2007t iBoccmet fe Barratll2007l) . The results are usually interpreted in terms 
of the Navier boundary condition, 



Au = u-U = b — 
on 



(1.1) 



where the fluid velocity u minus the surface velocity U is proportional to the shear strain 
rate du/dn via the slip length b. Flow past smooth hydrophilic surfaces has been shown 
to be consistent wi th the no-slip hypothesis, but b can reach tens of nanometres for hy- 
droph obic surfaces ( Vinogradova fc Yakubov|2003tlCottin-Bizonne et "aLl l2005l: iJolv et all 
2006| . Hydrophobicity can be significantly amplified by roughness and can reduce fr iction 



due to trapped nanobubbles (jVinogradova et aLlll995l ; ICottin-Bizonne et a/.ll2003r) . Ex- 



treme hydrophobicity can be gener ated with well-controlled textures riQuere!|2005l), lead- 
ing to a many-micron slip lengths (|Ou fc Rothsteinll2005t I Joseph et al\ 120061 : IChoi et al 
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2006) and very fast transport of water through microchannels. The strong anisotropy of 
such surfaces, however, can limit the validity of Eq. (|1.1|L 

The possibility of transverse flow over a grooved no-slip surface , perpend i cular to an 



ine possibility 01 transverse now over a grooved no-snp surrace , perpendicular to an 
applied shea r stress, has been analyzed by IStroo"ck~ et all (|2002al ). lAidaril J2002h. and 



Wand (l2003h and exploit ed for chaotic mixing in microfluidic devices by IStroock et al 



( 200261 ). In this context. IStroock et all (|2002al ) expressed the permeability k of a thin 
(parallel-plate) microchannel with one grooved and one flat surface in terms of an effective 
slip-length tensor, b = {fry}, defined by a generalized Navier boundary condition 



Au = u U = b (n • Vu) 



(1.2) 



and Stone et all ( 2004 ) expressed the velocity profile in terms of b. This elegant con- 
struction relating permeability to slip, however, assumes that the global flow has the 
same anisotropy as the grooved surface ( i.e. k and b are coaxial). This is generally no t 
the case with multiple textured surfaces ( Wang 2003 ). curved walls ( Einzel et al. 199Cll ). 
obstacles in the flow, etc., an d we are not aware of any other use of the tensorial relation 
(|1.2[) . Notably, |WangJ (|2003) considered flow between misaligned, grooved plates using 
(|1.2[) in component form but deemed the solution 'too tedious to reproduce'. We shall 
see that this problem and others have very simple solutions in tensorial form. 

In this article, we propose the use of (jl.2p as a local boundary condition for any surface 
whose texture perturbs fluid flow on length scales much smaller than the geometry. We 
begin in section [2] by discussing a general boundary condition relating slip velocity to 
normal traction via an interfacial mobility tensor. To illustrate its use, we derive exact 
solutions for two types of laminar flow between textured parallel plates (which can also be 
superimposed) : (i) shear flow due to moving plates in section [31 and (ii) pressure-driven 
flow in section [4] We close in section [5] by suggesting further applications. 



2. Theory 

2.1. The interfacial mobility tensor 

Although Equation (|1.1|) is the most common ly used bound ary condition for hydrody- 
namic slip, it is not widely appreciated that Navierl ( 18231 ) also postulated the more 
general relation, 

Au = Mt (2.1) 

where r is the local shear stress (normal traction) and M is a constant interfacial mobility 
(velocity per surface stress). For a Newtonian fluid, r = -qdu/dn, this reduces to (|1.1[) 
with b = Mr], where rj is the viscosity. Molecular dynamics simulations have shown that 
(|2.ip with constant M is more robust than (|1.1|) with constant b, since the fluctuating 
slip velocity correlates better with the shear stress (normal forces) than with velocity 
gradients very close to the surface ( Hess fe Loosdll989l : Bocquet fc Barratl 20071 ). 



A natural generalization of the slip condition (J2JJ) is 

Au = M(n-er) (2.2) 

where n ■ cr — f„ is the fluid force (normal traction) on the interface, er is the local 
stress tensor, and M is an interfacial mobility tensor. As shown in Figure [U the effec- 
tive slip vector is generally misaligned with the force vector for an anisotropic surface. 
Equation (| 1 . 2[) is recovered with b = M77 in the case of a Newtonian fluid of viscosity 77. 
For anisotropic surfaces, the mobility is a second-rank tensor M = {My}, whether the 
averaging of surface fo rces occurs over microstru ctural or molecular heterogeneity. 



As with scalar slip (jBocquet fc Barratl 12007). the tensorial slip boundary condition 
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Figure 1. Physical picture of tensorial slip. The normal traction f n exerted by the fluid on 
an anisotropic surface produces an effective slip velocity Au = Mf„ in a different direction. At 
the molecular level, the interfacial mobility tensor M can be related to trajectories of diffusing 
interfacial particles, such as the one shown. 



(|2.ip can be justified in various ways. At the microstructural level, grooved surfaces (with 
or without scalar slip) have effective tensorial slip coefficients, which can be explicitly 
calculated fo r simple geometries, if the groov es vary on much smaller scales than the 



fluid domain (jStroock et aLll2002atlWa ng 2003). At the molecular level, nanoscale surface 



anisotropy has a similar effect, but due to statistical interactions. 

A possible starting point for molecular modeling is a tensorial Einstein relation, D = 
MfcT/S 1 , relating the M to the 'interfacial diffusivity' per unit area S, by analogy with 
the theory of Brownian motion. This yields the statistical formula 



Mij = ^ Hm | Cov(Aa;,-(t), A Xj (t)) ( 2.3 ) 



where Ax(t) = x(t) — Ui is the fluctuating position of an interfacial fluid molecule, in a 
frame moving with the mean surface velocity (Fig. [T]), where the 'interface' may include 
molecules distinct from the bulk fluid, e.g. in a vapor phase. The idea of observing thermal 
diffusion near a surface to infer its slip length has also been exploited in experiments by 
Jolv et al. (2006). The mobility formula (|2.3|) can also be recast in a tensorial Green-Kubo 



M v = J7^ <aCov(«i(0), «,(«)) (2.4) 



n 



form, 

S_ 

kT , 

where v(i) = dAx/dt. These formal expressions assume convergence in the thermody- 
namic limit (taken before t — > oo) for molecular trajectories exploring the interfacial 
region on scales much larger than the surface heterogeneity. In that case, via the co- 
variance matrix, M is s ymmetric, positive definite and thus invertible. As noted by 



Bocquet fc Barrati (|2007l ). the inverse mobility, or friction tensor F = M , also has a 
tensorial Green-Kubo representation, as t he integral of the auto-co rrelation function for 
forces exerted by the fluid on the surface ( Bocquet fc Barrati 1994 ). 



2.2. General properties of 'M 

Regardless of its microscopic justification, we suggest adopting ()2.2|) as a general, interfa- 
cial constitutive relation for continuum mechanics. As with its bulk counterpart relating 
the stress and deformation rate, its form can be either derived from microscopic models 
or simply postulated and fit to experimental data, subject to certain constraints dis- 
cussed below. For a general 'nonlinear interface', the mobility tensor M could depend on 
the surface forces, as well as internal degrees of freedom, such as the local orientation of 
surface molecules or deformable microstructures; for example, hinge- like structures could 
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lead to different slip in opposite directions. For permeable surfaces with n ■ Au ^ 0, the 
mobility tensor may be represented by a 3 x 3 matrix with tangential-normal couplings, 
a possibility which has not been considered before to our knowledge. 

Here, we will focus on the simplest case of impermeable, macroscopically homogeneous, 
linear interfaces, where M is a constant 2x2 matrix in a local suitable coordinate system 
of the tangent plane. Below we will refer to the mobility tensor as defining the 'texture' 
of a surface up to a rotation, which sets the 'orientation'. The tensor formalism allows 
us to easily change the orientation of a texture, once a problem has been solved in terms 
of mobility tensors for a given geometry. The mobility simply transforms as 

MhS»MS_ 6 where Sg = ( °° S0 a (2.5) 

V -sinfl cos 9 J K ' 

is a matrix rotating the tangent plane by an angle 9. 

We also consider 'passive' surfaces, which do not transfer energy to the fluid. In that 
case, enforcing a positive rate of work wj on the slipping interface (M ^ 0), 

w/=f n -Au = f„-Mf n >0 (2.6) 

for any loading f n = n-cr implies that M must be positive definite. This argument is sim- 
ilar to the constra int of positive entropy pr oduction at a slipping boundary in irreversible 
thermodynamics ( Heidenreich et aZll2007t) . The statistical arguments above lead to the 



same conclusion, e.g. since the diffusivity D is positive definite for a passive surface. 
The eigenvectors of M correspond to special directions along which fluid forces do not 
produce transverse slip, and the (positive) eigenvalues are the corresponding directional 
mobilities. Since positive definite matrices are invertible, the boundary condition can also 
be expressed as f„ = F Au in terms of the (coaxial) friction tensor, F = M _1 . 

Diagonalization allows us to relate M to the position of the slip plane in (|2.2p . which is 
independent of the force f n . In each eigendirection ej, the tensorial boundary condition 
(|2.2p reduces to the scalar case (12.11). and the eigenvalue M i depends on the (arbitrary) 



choice of slip plane in the usual way ( Bocquet fc Barrati2007 ): for a Newtonian fluid (jl.ip . 



the slip-length hi — Miij is the position of the slip plane, relative to the (unique) depth 
of no slip extrapolated from a homogeneous bulk shear flow. By appropriately shifting 
the eigenvalues {Mi}, the same slip plane can be chosen for all directions. The mobility 
tensor is then constructed from the spectral decomposition, M = SMS -1 , where M is 
the diagonal matrix of eigenvalues and S the matrix of column eigenvectors. 



2.3. Symmetric mobility tensors 

Although we will make no further assumptions in our analysis below, a constant mobility 
tensor is usually symmetric, My = Mji, as in the statistical formulae above. This is 
also the case for the eff ective slip tensor der i ved by aver aging linear Stokes flows over 
grooved no-slip surfaces (jStroock et al. 2002 b; IWan g 2003). More generally, symmetry of 
M exemplifi es the widely used Onsager-Casimir re l ations of linear response near thermal 
equilibrium (jBocquet fc Barratlll994l lAidarilEool iHeidenreich et aZ.I[2007h . 

A 2 x 2 interfacial mobility matrix, which is symmetric and positive definite, has some 
useful mathematical properties. There always exists a rotation of the orthogonal (x, y) 
coordinate system of the tangent plane Sg , which diagonalizes the mobility, 



M : 



M|| 





M, 



_ f M || cos 2 9 + Mj_ sin 2 9 (M\\ - M± ) sin 9 cos 9 
(M|| - Mj_) sin 9 cos 9 M\\ sin 2 9 + M± cos 2 9 



(2.7) 
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Figure 2. Sketch of a fluid region \z\ < h/2 between upper (+) and lower (— ) parallel plates 
with arbitrary textures (grooves, surface coatings, bubbles, etc.). The texture length scales are 
much less than the gap h, so each surface has a well defined slip-length tensor, b + and b~, with 
eigenvalues 6^ (and irjj) in the fastest (and slowest) slipping directions indicated. The plates 
move at relative velocity U and/or a uniform pressure gradient — Vp is applied. In addition to 
the usual no-slip parabolic Poiseuille flow P and linear shear flow S, there are superimposed 
slip-driven plug-flow A and shear-flow B in different directions. 

where M\\ > M± > are the eigenvalues. The decomposition M = M||e||eT + M±e±e^_ 
shows that Au is a linear superposition of scalar slip in the eigendirections. 

Regardless of the complexity of the texture, as long as Equation (|2.2[) holds at the ge- 
ometrical scale with a symmetric, positive definite M, there exist orthogonal directions 
on the surface, en = Sgx and e± = Sgy, along which there are no transverse hydrody- 
namic couplings, e± ■ Men = 0. The mobility for 'forward' slip aligned with forcing in a 
particular direction e — S^x is given by e • Me = My cos 2 (8 — (f) + M± sin 2 (9 — <j>) and 
is bounded by the eigenvalues, M± < e ■ Me < Mm . The 'fast' axis of greatest forward 
slip (6 = 0) is always perpendicular to the 'slow' axis of least forward slip (6 = ir/2). 

2.4. Simple examples 

In the following sections, we focus on passive linear interfaces and Newtonian fluids, 
described by the Navier-Stokes equations 



In that case, all the properties of M above are inherited by the slip-length tensor, 
b = Mr;, with eigenvalues, 6|| = My 77 and 6^ = M±_rj. We also assume impermeable, 
macroscopically homogeneous surfaces, for which b is a constant 2x2 matrix. 

To illustrate the use of slip tensors, we consider the geometry in Fig. [2] where the fluid 
is confined between flat plates at z = ±h/2 moving at velocities XJ^ (this section) or 
forced by a pressure gradient (next section). Each plate has a fine texture (varying on 
scales h) and exhibits uniform tensorial slip, 




and V ■ u = 



(2.8) 



u = U± T b±^ 

oz 



for z = ± — 
2 



(2.9) 



where the slip-length tensors, b and b , are represented by constant, positive definite 
(but not necessarily symmetric) 2x2 matrices in the (x, y) coordinate system. 



6 



M. Z. Bazant and O. I. Vinogradova 



3. Example: Shear flow 

3.1. General solution 

The simplest solution of (|2.8p - (|2.9p corresponds to laminar shear flow between two moving 
textured plates, shown in Fig. [21 In terms of the depth-averaged velocity U = (U + + 
U~)/2 and relative plate velocity v = U + — U~, we can express the solution as 



u = U 



(B. + I)f 



- 2 (3-D 



where A s and B s are dimensionless 2x2 matrices with the following physical interpreta- 
tions. The first term in (|3 . 1 [) describes a slip-driven plug flow in the A s v direction, and 
the second describes a slip-driven linear shear flow in the B s v direction. Substituting 
([3~Tj) into (HU), we find 

A^-Dfl + C)- 1 and B = (I + C)- 1 -I, (3.2) 

where 

hC = h + +\>- and /iD = b+ b . (3.3) 

The slip-driven plug flow vanishes (A s = 0) only if the textures are the same (D = 0), 
and slip-driven shear flow always occurs (B s ^ 0, if b + ^ or b~ ^ and thus C / 0). 

The solution (|3.1[) - (|3.2p exists for any b ± , as long as I + C is invertible; this is ensured 
for passive surfaces, since ^ and C are positive definite, and possible for some active 
surfaces. In the typical case of symmetric b ± , the solution can be expressed in terms 
of the texture orientation angles 6 I± and slip-length eigenvalues, — M^rj and b± — 

M±i~i using (|2.7j) . This can be easily accomplished in the following general situations by 
diagonalizing A s and B s . 

3.2. Aligned but different textures 

We first consider 'aligned' surfaces with the same orientation 0^ = 8, but arbitrary 
slip-length eigenvalues: 



b 



The coefficient tensors (|3.2[) are then diagonalized by the same rotation matrix 
and the eigenvalues 

result from scalar slip in the eigendirections. 

3.3. Identical but misaligned textures 
Next we consider identical textures with arbitrary orientations, 6' ± = 9 ± AO: 

b± = S e± ( ^ £ \ S_ e± (3.8) 
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By expressing the sum and difference matrices ([3T3J) as 

hC = 28, ( & " C ° s2 M + ^ Sin2 A0 fc . 2 ° fc 2 Afl ) S s (3.9) 
\ 6|| sm 2 A9 + b± cos 2 A0 ' 



/ (&||-6 ± )sin2A0 

1 b " 1 (6|| -6 ± )sin2A0 



we find A, = and 



bii cos 2 A6»+6j_ sin 2 A8 



(3.10) 



B-=S„| ^ /2+fc||C °7^ Si " A9 b „sin 2 A^cos 2 A |S-,- CUI. 

U /i/2+&|| sin 2 A9+6j_ cos 2 Ad 

The slip-driven plug flow vanishes by symmetry , and the slip-driven shear flow coefficient 
B s is diagonalized by S§, where 9 is the angle that bisects the orientation angles. As 
expected by symmetry, if the two textures are the same (but misaligned) , then shearing 
in this direction cannot produce any transverse flow. 



4. Example: Pressure-driven flow 

4.1. General solution 

Another simple solution to (|2.8[) - (|2.9[) describes steady, laminar flow in response to an 
applied pressure gradient, g = — Vp = g x x + g y y, between stationary textured plates in 
Fig. O We express the solution in the form 




>-'! 



A,. 134 g (4J) 



where A p and B p are dimensionless 2x2 matrices. In spite of surface anisotropy, the 
velocity is horizontal [z ■ u = 0) and varies only in the vertical z direction, due to 
translational invariance. The solution (|4.1|) is a linear superposition of three terms: The 
first is the familiar parabolic profile of Poiseuille flow in the g direction between parallel 
no-slip planes; the second is a slip-driven plug flow in the A p g direction; the third is a 
linear shear flow in the B p g direction, which arises only if b + ^ b~. Substituting (|4.1[) 
into (HH) , we find 

A 1 , = C-D(I + C) _1 D and B p =(I + C)- 1 D (4.2) 



in terms of the sum and difference tensors defined in 

In the limit of no slip on the upper surface b + = 0, our solution reduces to that 
of Stone et al. I (l2004 . In that case, hC = —KD = b , the coefficient tensors, A p and 



Bp, and the permeability K are all coaxial with the slip- length tensor b~ of the lower 
surface. Here, we analyze more general situations where the upper and lower surfaces 
have different slip tensors. 

For symmetric h^, we can diagonalize A p and B p in the same simple situations con- 
sidered above for shear flow. In the case of aligned but different slip tensors (|3.4|) , the 
coefficient tensors (|4.2[) are diagonalized by the same rotation matrix: 
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and the eigenvalues 

a /l+ l \ b+ + b- +4h- 1 b+b- _ ,_. 

A s (6+,6 )= , + b+ + b - -d fl.(6+6 ) = - - fo+ - fc _ (4.5) 

result from scalar slip in the eigendirections. There are several simple cases: (i) If the 
surfaces are isotropic, b = 6 ± I, then A p = A S I and B p = B s l; (ii) if the surfaces have 
the same slip tensors, b + = b = b, then A p = 2h~ 1 h and B p = 0; (hi) If the upper 
surface has no slip, b + = and b = b, then A s = —B s = bj(h + b) , or more compa ctly 
A p = —Bp = b(hl + b) _1 , which reduces our solution to that of lStone et al. I (|2004h for 
one textured surface. 



For identical but misaligned textures (|3.8p , we find A p = C and B p = 0, using 
and (|3.10[) . Now the slip-driven shear flow vanishes by symmetry. The slip-driven plug 
flow is proportional to the average slip-length tensor and diagonalized by S§, where 9 is 
the angle that bisects the surface orientation angles. As expected, a pressure gradient in 
this direction cannot produce any transverse flow, if the two textures arc the same. 

4.2. Permeability 

In many situations, one is more interested in the depth-integrated total flow rate in a 
given direction, rather than the velocity profile. In linear response, the depth-averaged 
velocity u is proportional to the applied pressure gradient, 



1 



h/2 



u= — / udz = Kg. (4-6) 
h J -h/2 

via the permeability tensor n. For the anisotropic Poiscuillc flow (|4.ip . this integral is 
easily performed to obtain 

h 2 

K = j^K, where K = I + 3A p (4.7) 

is the dimensionless permeability, scaled to its value without slip. The permeability is 
generally enhanced by slip-driven plug flow in the direction A p g. (The slip-driven shear 
flow does not affect the permeability, although it contributes to mixing and dispersion.) 

The results above for A p in various special cases can be extended to K 1 since the two 
tensors are coaxial: 

K = I + 3Ap = S 0K ( ^11 ^ \ S„ 9k (4.8) 

where Sg K diagonalizes A p and K. For aligned but different textures (13. 4p . the perme- 
ability clearly has the same orientation as the textures, 9k — 9, and its eigenvalues, 
if II = K s (bt , fejp and K± = K s {b\ 1 &][) correspond to analogous cases of scalar slip, 

For identical but misaligned textures p.8p . the permeability is orientated with the mean 
angle 9k = 9 — (9 + + 9~)/2 with eigenvalues given by 

6 
h 



if II = 1 + ^ (6 !t cos 2 AO + b ± sin 2 AO) (4.10) 
6 



K ± = 1 + ^ (6|| sin 2 AO + b ± cos 2 AO) (4.11) 

where A9 = (9 + — 9~)/2. If A9 — 7r/4, then the permeability is isotropic, K = KI with 
K = K\\ =K ± = l + (3/h)(b\\+b±). 
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Microfluidic devices often contain thin channels of rectangular cross section with par- 
allel side walls at y = ±L with L 3> h. In that case, the mean downstream permeability 
of the channel, k x = (h 2 /12i])K x , defined by Ux — ^xQx i CcLH be easily derived from 
the permeability tensor k defined in Eq. (|4.6p . Ignoring departures from Poiseuille flow 
within 0(h) of the side walls, the constraint of vanishing transverse flow, u y — 0, is 
maintained by an induced transverse pressure gradient, g y — — (n yx /K y y)g x , which drives 
an additional anisotropic Poiseuille flow. Superimposing these flows, we obtain 

K. = K„-!^ = *>M=M± (4.12) 
K K K 

I\yy l\yy I\yy 

The channel permeability can also be interpretted in terms of an effective downstream 
slip length b x defined by K x = 1 + (6/h)b x , although this obscures the true tensorial 
nature of the hydrodynamic slip. 



5. Conclusion 

Our solutions for anisotropic flows between textured plates may be useful in interpret- 
ting experiments and simulations. As in the case of scalar Poiseille flow, bulk velocity 
profiles can be fitted to the theory to systematically extract boundary effects of slip- 
page and assess the validity of the tensorial slip hypothesis. Our results also allow the 
local slip tensors to be determined by global measurements, such as the permeability of 
a textured channel or the force required to shear textured plates, as a function of the 
surface orientations. In such measurements, departures from our predictions could be 
used to isolate nonlinear, inhomogeneous, or non-symm etric slip response, e.g. due t o 
nanobubble defor mation at superhyd rophobic surface ( Sbragaglia fc Prosperettil 2007 ). 
surface curvature ( Vinogradova 1995t) . or variable channel width ( Lauga et al. f l2004l) . 

More generally, our calculations illustrate the power of the tensor formalism to cap- 
ture complicated effects of textured surfaces, while preserving simple fluid domains. The 
general boundary condition (12. 2p may be useful for analytical or numerical calculations 
in many other situations, such as lubrication flows between text ured gears, spreading o r 
drainage of thin films, dispersion and mixing in groo ved channels (ptroock et 
sedim ention of text ured particles (jLecoa et and electrokinetics of patterned sur- 

faces (|Aidarill2002f) . Transverse spatial couplings could also be added to existing tensorial 
(but isotropic) slip boundary conditions for fluids with internal degrees of freedom, such 
as liquid crystals and polymer melts; anisot ropic surface texture ca n influence molecu- 
lar orientations and thus the effective slip (jHeidenreich et all 120071 ). which could have 
interesting consequences for theory and applications. 

The authors gratefully acknowledge the hospitality of ESPCI and support by the Paris- 
Sciences Chair (MZB) and Joliot Chair (OIV). 
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